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Abstract: Drought acutely affects economic sectors, natural habitats and communities. Understanding the 
past spatial and temporal patterns of drought is crucial because it facilitates the forecasting of future 
drought occurrences and informs decision-making processes for possible adaptive measures. This is 
especially important in view of a changing climate. This study employed the World Meteorological 
Organization (WMO)-recommended standardized precipitation index (SPI) to investigate the spatial and 
temporal patterns of drought in Zambia from 1960 to 2016. The relationship between the occurrence of 
consecutive dry days (CDD; consecutive days with less than 1 mm of precipitation) and SPI was also 
investigated. Horizontal wind vectors at 850 hPa during the core of the rainy season (December—February) 
were examined to ascertain the patterns of flow during years of extreme and severe drought; and these 
were contrasted with the patterns of flow in 2007, which was a generally wet year. Pressure vertical 
velocity was also investigated. Based on the gamma distribution, SPI successfully categorized extremely 
dry (with a SPI value less than or equal to —2.0) years over Zambia as 1992 and 2015, a severely dry (-1.9 
to —1.5) year as 1995, moderately dry (—1.4 to —1.0) years as 1972, 1980, 1987, 1999 and 2005, and 26 near 
normal years (0.9 to 0.9). The occurrence of CDD was found to be strongly negatively correlated with 
SPI with a coefficient of —0.6. Further results suggest that, during wet years, Zambia is influenced by a 
clockwise circulating low-pressure zone over the south-eastern Angola, a second such zone over the 
northern and eastern parts, and a third over the Indian Ocean. In stark contrast, years of drought were 
characterized by an anti-clockwise circulating high-pressure zone over the south-western parts of Zambia, 
constraining precipitation activities over the country. Further, wet years were characterized by negative 
pressure vertical velocity anomalies, signifying ascending motion; while drought years were dominated by 
positive anomalies, signifying descending motion, which suppresses precipitation. These patterns can be 
used to forecast drought over Zambia and aid in strategic planning to limit the potential damage of 
drought. 
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1 Introduction 


The World Meteorological Organization (WMO) defines a meteorological drought as an extended 
period of dryness in the natural climate cycle compared to the 30-year average of the region under 
study (Monacelli, 2005). Exploring patterns of drought is of interest because of its detrimental 
effects on socio-economic activities. These impacts have been well documented for different parts 
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of the world in recent years (Santos et al., 2010; Udmale et al., 2014). Over Africa, a recent study 
conducted by Uhe et al. (2017) showed that Kenya suffered severe drought following low 
precipitation coupled with high temperature towards the end of 2016. 

In the face of a changing climate (Field et al., 2012), the intensity and frequency of drought are 
projected to increase with an increase in radiative forcing, especially by the end of the 21" 
century (Taylor et al., 2012; Zhao and Dai, 2015). These changes will lead to damaging impacts 
on natural ecosystems, cause dire shortages of basic needs such as water and food and lead to an 
increase in heat waves, resulting in wild fires (Dai, 2011). Statistical links between climate 
change and occurrences of drought have been documented (Burke, 2011). Although it has been 
noted that below normal precipitation is the main driver of drought (Trenberth and Fasullo, 2009; 
Zhao and Dai, 2015), a warming climate increases evapotranspiration, which may lead to more 
pronounced droughts (Scheff and Frierson, 2012). Southern Africa is particularly vulnerable to 
climate-change-induced drought because of its dependence on climate-driven sectors such as 
rain-fed agriculture, and services provided by ecosystems (Loomis, 1999). 

As part of southern Africa, Zambia has been experiencing the impacts of drought for some time 
now. For example, Hill (2016) attributed the drying up of Lake Kariba to drought. Lake Kariba, the 
world's largest man-made lake by volume, contributes significantly to hydroelectricity, fisheries, 
and agricultural activities (Muchuru et al., 2014). It's drying up has therefore been a source of 
concern to policy makers and has attracted regional and international attention (Hill, 2016; New 
Yorker, 2016). Invariably, during years that Zambia has experienced drought, Lake Kariba has 
experienced dangerously low water levels, leading to a reduction in power generation and 
consequently adversely affecting socio-economic activities in both Zambia and Zimbabwe (Hill, 
2016). Although there is some documentation on the impacts of drought in Zambia, patterns of 
drought occurrences, and circulation mechanisms driving the drought remain poorly understood. 

Given the negative effects of drought, it follows that the primary protection against this 
phenomenon is a proper understanding of its past patterns in order to understand its probable 
future occurrence(s). In examining past drought patterns, diverse methodological approaches have 
been used around the world. Some studies (e.g., Okoola, 1999; Ogwang et al., 2012) have utilized 
composite analysis to examine atmospheric mechanisms affecting precipitation during drought 
years. 

In this study, the standardized precipitation index (SPI) was used to explore the patterns of 
drought because of its convenience in categorizing them (McKee et al., 1993; Hayes et al., 2011; 
Ngarukiyimana et al., 2017). SPI is based on precipitation data only, and it runs a comparative 
analysis of annual average precipitation and transforms the data to a normal distribution in order 
to avoid errors. The transformation to a normal distribution is performed using an 
equal-probability transformation after fitting the precipitation data to a gamma distribution 
(Zargar et al., 2011). The average is then fitted to zero with above-zero values denoting relatively 
high precipitation and below-zero values denoting drought. 

The main goal of this study is to investigate spatial and temporal patterns of drought in Zambia 
over a 56-year period (1960-2016). As in Zargar et al. (2011), drought was used in this study to 
represent a stochastic natural hazard, which is triggered by rainfall failure associated with 
persistent synoptic mechanisms, affecting water resources, agriculture, and socio-economic 
activities. As well as precipitation itself, circulation mechanisms related to extremely and severely 
dry years are also investigated in this study. The findings highlight mechanisms that must 
carefully be monitored in the future for strategic planning to avert losses associated with drought. 


2 Materials and methods 


2.1 Study area 


Zambia is located in the southern parts of Africa (longitudes 21°30’E to 34°00’E and latitudes 
18°00’S to 08°00’S). Generally, much of the landmass in Zambia lies between 910 and 1200 m 
above sea level (a.s.l.). Zambia can therefore be classified as a plateau with peak points (more 
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than 1800 m a.s.1.) being found along the Muchinga Escarpment (northeast) and the Mafinga hills 
(east) and lower points (less than 610 m a.s.1.) in the valleys of the Luangwa and Zambezi rivers. 

Zambia experiences a bimodal precipitation pattern with a distinct dry season from April to 
September and a wet season from late October to early April next year. Many studies (e.g., Tiffen 
and Mulele, 1994; Hachigonta and Reason, 2006; Libanda et al., 2015) have found statistical links 
between the onset/cessation of rainfall and the movements of the Inter Tropical Convergence 
Zone (ITCZ). The ITCZ is a low-pressure system that occurs when the north-eastern trade winds, 
moist air from the Congo basin, and the meridional south-eastern trade winds converge (Libanda 
et al., 2017a). This low-pressure zone triggers precipitation events. The ITCZ begins to affect 
Zambia from the north and north-west around October or November, causing the northern half of 
the country to experience increased precipitation while most parts of the southern half remain 
relatively dry. From December to February, the ITCZ travels southwards bringing precipitation to 
the southern half of the country before withdrawing to give a second precipitation peak in the 
northern half in March or April (Libanda et al., 2017b). Total annual average precipitation in the 
northern half of the country is 1400 mm while the southern half receives about 600 mm (Libanda 
et al., 2015). 


2.2 Data sources 


The datasets used in this study are summarized in Table 1. Specifically, Global Precipitation 
Climatology Centre (GPCC) precipitation data were used to calculate the occurrences of drought. 
This dataset is monthly, and it is gridded with a spatial resolution of 0.5°x0.5°, 1.0°1.0°, and 
2.5°x2.5° latitude by longitude. It is described in detail by Schneider et al. (2016) and is available 
through the Physical Science Division of the Earth System Research Laboratory 
(https://www.esrl.noaa.gov/psd/data/gridded). 


Table 1 Summary of datasets the used in this study 


Dataset Source Resolution Period used 


Physical Science Division, Earth System 
Research Laboratory 

ARC2 Climate Prediction Center 0.1°x0.1° 1983-2016 

2007 (wet); 1995 (severely 


GPCC 0.5°x0.5° 1960-2016 


Relative humidity, zonal European Centre for Medium-Range Weather 0.75°x0.75° dry); 1992 and 2015 (extremely 
and meridional winds Forecasts dry) 
Pressire verticalvelocit National Centers for Environmental 2007 (wet); 1995 (severely 
3 : ey Prediction/National Center for Atmospheric DIRLS dry); 1992 and 2015 (extremely 
(omega) 
Research dry) 
Topographical data Joint Institute for the Study of the Atmosphere 1°x1° 
and Ocean 


In this study, the maximum number of consecutive days with less than 1 mm of precipitation 
was investigated because it serves as a reliable measure of the evolution of seasonal drought 
(Duan et al., 2016). Originally, the theory of consecutive dry days (CDD) was proposed by Frich 
et al. (2002) in their study on variations and changes in extreme climate towards the end of the 
20" century. Subsequently, during the CLIVAR/GCOS/WMO workshop on indices and indicators 
for climate extremes, the Expert Team on Climate Change Detection and Indices (ETCCDI) 
adopted this index and included it in the 27 core climate change indices, which have since been 
employed by many studies quantifying climate change worldwide (Karl et al., 1999; Ongoma et 
al., 2017; Libanda and Ngonga, 2018). 

To investigate the evolution of CDD, we obtained daily rainfall data (Africa Rainfall Climatology 
Version 2; ARC2) from the archives of the Climate Prediction Center and the data are freely 
available at https://iridl.ldeo.columbia.edu. This dataset is a blend of rain gauge and satellite 
information. It is of high quality and is gridded at a spatial resolution of 0.1°x0.1° (Novella and 
Thiaw, 2013). The period covered for the calculation of CDD in this study is 1983 to 2016. We 
chose this period based on data availability. Unlike GPCC monthly data that covers a longer period 
of time, ARC2 daily rainfall is only available for the period 1983-—present. 
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Relative humidity (RH), and zonal and meridional winds were studied, on the basis of 
ERA-interim reanalysis data (https://www.ecmwf.int; Dee et al., 2011). These datasets are gridded 
at a resolution of 0.75°x0.75°. They are of high quality, and many studies (e.g., Nicholson and Grist, 
2003; Bordi et al., 2005; Libanda et al., 2016) have used them in different parts of the world. 

The topographical dataset used in this study is a product of the Joint Institute for the Study of the 
Atmosphere and Ocean (JISAO; http://jisao.washington.edu). It is a 1°x1° latitude-longitude 
resolution elevation dataset. 


2.3 Data analysis 
Using a gamma distribution probability density function, SPI was computed as: 


tear 
80) = Farr (>), (1) 


x 


where J (a) represents the gamma function; x is the amount of precipitation in mm (x>0); a is the 
shape parameter; and J is the scale parameter (6>0). Many studies (e.g., Dogan et al., 2012; 
Salehnia et al., 2017) have used this method for drought characterization in different parts of the 
world. 

The classification of drought years in this study was performed following McKee et al. (1993) 
and it is given in Table 2. The SPI has been calculated on an annual basis. 


Table 2 Classification of drought on the SPI scale 


Description SPI magnitude 
Extremely dry <2 

Severely dry —1.9 to -1.5 
Moderately dry —1.4 to-1.0 
Near Normal —0.9 to 0.9 


Drought is a duration phenomenon in that it refers to a time span, therefore, understanding the 
evolution of CDD is critical to better inferring the variability of drought. The occurrence of CDD 
and its relationship to seasonal drought was therefore explored using the Taylor diagram (Taylor, 
2001). This method is a convenient way of summarizing the statistical association between two or 
more variables, and has been used in many climate studies (e.g., Fotso-Nguemo et al., 2018; 
Ongoma et al., 2018). The parameters involved are correlation, standard deviation and root mean 
square error. 

Composite analysis was employed to investigate atmospheric patterns influencing precipitation 
during drought years. Composite analysis is a useful technique that involves the collection and 
averaging of meteorological variables of interest. These composites are then used to interpret 
patterns associated with given scenarios (Okoola, 1999; Ogwang et al., 2012). Here, composites 
of wind flow at 850 hPa, pressure vertical velocity (omega) and RH were constructed. The years 
considered were extreme (1992 and 2015) and severe (1995) drought years. These were compared 
with composites of the most recent wet year, 2007. The aim of this was to investigate circulation 
anomalies associated with drought years in contrast to those associated with wet years. 

Following Singh and Rathor (1974), pressure was used in this study in place of altitude to 
investigate vertical motion. Mathematically, vertical velocity is given in pressure coordinates as: 


Dp 

á Dt’ @) 
where œ represents the Lagrangian rate of change in pressure with time (Pa/s); and Dp and Dt 
represent the time rate of change of pressure (Pa/s). Other studies (e.g., Libanda et al., 2016) refer 

the omega equation, a partial differential equation for vertical velocity given by: 

2 
Waor = ¢ by xVulE, “Avil xV al (3) 
op op op 
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The following definitions are used in Equation 3 (see Holton, 1992). The static stability is 
given by ø (rad/s, which is usually abbreviated as s"); y, is the horizontal del operator; y} is 


the horizontal Laplacian operator; f is the Coriolis parameter, with the unit of rad/s 
(approximately 10 */s, varying with latitude); vg is the geostrophic velocity vector (m/s); & is the 
geostrophic relative vorticity and is generally a scalar rotation thus dimensionless; ¢ is the 
geopotential (m); and p is defined as pressure (Pa). 

The main tools used for data analysis and presentation of results in this study are the R 
programming language, Climate Data Operators (CDO) and Grids Analysis and Display System 
(GrADS). The MDM (Meteorological Drought Monitoring) software was used to calculate SPI. A 
detailed description of this software is available in the work of Salehnia et al. (2017). 


3 Results 


3.1 Seasonal cycle of precipitation 


In this study, the core of the rainy season was observed (Fig. 1) to occur during the 
December—February (DJF) months with over 60% of cumulative precipitation received during 
these months. This observation was also made by Hachigonta and Reason (2006) who studied the 
inter-annual variability of dry and wet events over Zambia. 
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Fig. 1 Mean cycle of precipitation over Zambia for the period 1960-2016 based on GPCC (Global Precipitation 
Climatology Centre) data. The values were averaged over longitudes 21°30'E-34°00'E and latitudes 
18°00'S—08°00'S. 


3.2 Standardized precipitation index (SPI) 


Using the methods discussed in Section 2.2 of this paper, SPI results shown in Figure 2 indicate 
that Zambia experienced more dry years (28) than wet years (27) during the period 1960-2016. 
The most recent extremely dry event occurred in 2015. With reference to Table 2, the SPI 
categorized extremely dry (with a SPI value <-2.0) years over Zambia as 1992 and 2015, one 
severely dry (—1.9 to —1.5) year as 1995, moderately dry (—1.4 to —1.0) years as 1972, 1980, 1987, 
1999 and 2005, and 26 near normal years (—0.9 to 0.9). 

The monthly cycles of precipitation for extremely dry (1992 and 2015) and severely dry (1995) 
years were investigated to ascertain how they differ from a wet year (2007). The results are 
presented in Figure 3, which indicate that while precipitation started in October during both the 
extremely and/or severely dry years and the wet year, the overall magnitude was different, with 
less than 200 mm being received during the dry years and over 300 mm during the wet year, by 
the end of December. Similarly, the end of the rainy season was in late April or early May for both 
the dry years and wet year. However, the wet year was associated with more intense precipitation 
than the dry years. 
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Fig. 2 Standardized precipitation index (SPI) results over Zambia for the period 1960-2016 based on GPCC 
data. The values were averaged over longitudes 21°30'E—34°00'E and latitudes 18°00'S—08°00'S. 
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Fig. 3 Monthly cycles of precipitation for extremely dry years (1992 and 2015), severely dry year (1995) and 
wet year (2007) 


The observed extremely (1992 and 2015) and severely (1995) dry years were investigated in 
terms of their spatial distribution of precipitation during the core of the rainy season (i.e., DJF). 
We subtracted drought years from the long term mean (LTM) to clearly capture the departure of 
precipitation from the LTM (Fig. 4). Results indicate a wetter northern half of Zambia and a drier 
southern half. Generally, the DJF LTM precipitation ranged from 800 to 850 mm over much of 
northern Zambia and was less than 800 mm over the Livingstone region in southern province. 
During the severely dry year 1995, negative departure from the LTM was experienced over much 
of the country within the range of 0-350 mm. However, positive anomalies of approximately 100 
mm were observed around Mbala and the Lake Mweru Wantipa region (Fig. 4b). Figures 4c and d 
show the spatial variability during the extremely dry years of 1992 and 2015. The negative 
anomalies varied between 0 and 500 mm over much of Zambia. It is interesting to note however, 
that, similarly to DJF in 1995, positive anomalies of approximately 100 mm were observed in 
Lake Mweru and Lake Mweru Wantipa region of northern Zambia during the DJF in 2015. 


3.3 Evolution of consecutive dry days (CDD) 


In this study, the occurrence of CDD was investigated in terms of their statistical association with 
SPI (Fig. 5). To ease the computation, we removed the influences of dispersion by standardizing 
both CDD and SPI. Results shown in Figure 5 indicate a correlation coefficient of 0.6. The 


201903.00241v1 


chinaXiv 


ChinaXivA FEAF 


JOURNAL OF ARID LAND 


23°E 27°E 31°E 23°E 27°E 31°E 


ee 10°S 10°s 10°S 


Precipitation 
anomaly (mm) 
14°S 14°S 14°S 


18°S = 
23°E 27°E 
23°E 27°E 31°E 


18°S 
10°S 10°S $ _4 10°S 
14°S 


14°S 14°S 


18°S 18°S 
23°E 27T°E 31°E 


18°S 


23°E 27T°E 31°E 


Fig. 4 DJF (December—February, the core of the rainy season) precipitation anomalies over Zambia, based on 
GPCC data. (a), long term mean (the period 1960-2016); (b), severely dry year 1995; (c), extremely dry year 
2015; and (d), extremely dry year 1992. 
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Fig. 5 Taylor diagram of consecutive dry days (CDD) versus SPI for the period 1983-2016 based on Africa 
Climatology Version 2 daily rainfall data. The values were averaged over longitudes 21°30'E-34°00’E and 
latitudes 18°00’S—08°00’S. RMS, root-mean-square. The Taylor diagram has three axes: correlation coefficient 
(azimuthal position; black dotted line), standard deviation (radical distance from the origin; red dotted line) and 
centered RMS difference (distance from the reference point (green square); green dotted line). 


centered root-mean-square (RMS) difference between CDD and SPI was given relative to the 
distance of the green square on the x-axis. The green contours show the RMS values. In this case, 
the RSM is 1.76. The standard deviation of CDD is proportional to the radial distance from the 
origin and is here equal to 0.84. These results are in line with theoretical climatology, i.e., the 
more CDD are experienced in a season, the higher are the chances of experiencing drought (Duan 
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et al., 2016). 

A time series of CDD and SPI is shown in Figure 6. In line with the results from the Taylor 
diagram analysis, these results show that the occurrence of CDD was negatively correlated with 
SPI. Although the difference between CDD and SPI is minimal in some years (e.g., 1988), it is 
notable that there is a larger difference during extremely dry (i.e., 1992 and 2015) and severely 


dry (1995) years. 


1983 1986 1989 1992 1995 1998 2001 2004 2007 2010 2013 2016 


Fig. 6 Temporal variability of CDD and SPI for the period 1983—2016 based on Africa Climatology Version 2 
daily rainfall data. The values were averaged over longitudes 21°30'E-34°00'E and latitudes 18°00'S—08°00'S. 


3.4 Circulation patterns 


Anomalies in horizontal wind vectors at low-pressure levels are a key ‘triggering’ mechanism for 
precipitation. Convergence triggers the ascent of air and consequently, condensing of air. In 
contrast, divergence favors subsidence of air and consequently, drying of air (Ogwang et al., 
2012). In this study, the horizontal wind vectors at 850 hPa during the core of the rainy season 
(i.e., DJF) were investigated to ascertain the patterns of flow during years of extreme and (1992 
and 2015) and severe (1995) drought. These were then contrasted with the patterns of flow in 
2007, a generally wet year over Zambia. Results suggest that, during wet year, Zambia is 
influenced by a clockwise circulating low-pressure zone over south-eastern Angola, a second such 
zone over northern and eastern provinces, and a third over the Indian Ocean. These low-pressure 
zones are indicated with a red-colored 'L' in Figure 7a. In stark contrast, years of dryness (Figs. 7b 
and c) are characterized by an anti-clockwise circulating high-pressure zone over the 
south-western parts of Zambia. These high-pressure zones are marked with a red-colored 'H' on 
Figures 7b and c. Figure 7c was generated by averaging the DJF values for 1992 and 2015. 
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Fig. 7 Horizontal wind vectors at 850 hPa during the core of the rainy season (i.e., DJF) for the wet year 2007 
(a), the severely dry year 1998 (b), and the extremely dry years 1992 and 2015 (c), based on ERA-interim 
reanalysis data. The letters 'L' indicates the low-pressure zone while 'H' indicates the high-pressure zone. 


Mean meridional cross section of DJF pressure vertical velocity (omega; Fig. 8) was also 
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investigated with the aim of understanding the pressure vertical velocity associated with 
extremely dry years (1992 and 2015) and severely dry year (1992). Again, they were contrasted 
with the wet year 2007. Results show that 2007 (Fig. 8a) was characterized by negative anomalies 
which indicate rising motion especially over parts of north-western and Copperbelt provinces, i.e., 
longitudes 21°30’E—28°00'E. In contrast, dry years (Figs. 8b and c) were dominated by positive 
anomalies, which indicate sinking motion and suppression of precipitation. 
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Fig. 8 Mean meridional cross section of DJF pressure vertical velocity (omega) anomaly (x10~ Pa/s) for the 
wet year 2007 (a), the severely dry year 1995 (b), and the extremely dry years 1992 and 2015 (c), based on 
ERA-interim reanalysis data. The red upward arrow in (a) indicates rising motion while the downward arrows in 
(b) and (c) indicate sinking motion. 


To further understand circulation patterns associated with drought years, we studied the amount 
of water vapor in the atmosphere (Fig. 9). This is usually expressed as a percentage (%) as RH 
(Grinn-Gofron and Bosiackaet, 2015). Previous studies have found strong statistical links 
between the occurrence of rainfall and high RH over a region. For instance, in a study of the 
2010-2011 drought over Kenya, Ongoma et al. (2015) found positive anomalies of RH over the 
western part of the study area being positively correlated with above normal rainfall. Similarly, in 
this study, the spatial patterns of RH were explored. Results indicate that the wet year 2007 was 
associated with high RH values especially over north-western and Luapula provinces (Fig. 9a). 
Notably, the RH exhibited a reduction from the north to the south of the country and was 
positively correlated with rainfall patterns (Fig. 4a). During severely dry year (1995), most parts 
of the country experienced a reduction in RH with only north-western, Copperbelt and Luapula 
provinces experiencing an average of 60%. The RH was less than 60% in the rest of the country. 
During extremely dry years (1992 and 2015), the whole country experienced a maximum average 
RH of only 50%. 


4 Discussion 


In this study, we found that while the observed occurrences of drought (Fig. 2) seem sporadic, 
Zambia experienced a 5-year period of drought from 1990 to 1995, including the extremely dry 
year of 1992 that led to critical food shortages coupled with myriad societal problems (Tiffen and 
Mulele, 1994). It is notable that more precipitation was observed over Luapula, northern and 
north-western provinces of Zambia, with southern and parts of western provinces receiving the 
lowest amounts (<200 mm). This observation was also made by Libanda et al. (2015) and by 
Hachigonta and Reason (2006) who described it as 'a downward gradient of precipitation from the 
north to the south of the country’. This downward gradient can be attributed to the influence of a 
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Fig. 9 Mean DJF relative humidity (RH) for the wet year 2007 (a), the severely dry year 1995 (b) and the 
extremely dry years 1992 and 2015 (c) 


persistent low-pressure zone over eastern Angola (Fig. 7a), which distributes moisture to the 
northern half of Zambia. It is important to note that Lake Kariba, Zambia's main source of 
hydroelectricity, is located south of Lusaka, a region exhibiting stronger drought patterns. Zambia 
sources more than 95% of its electrical -energy needs from hydroelectricity. In recent years, 
whenever the country has experienced drought, Lake Kariba has had dangerously low water 
levels, leading to critical reductions in power generation and, as a consequence, adversely 
affecting all socio-economic activities in both Zambia and neighboring Zimbabwe (Hill, 2016). 
These results can therefore inform decision-making processes in water resource management. 

The extremely dry (1992 and 2015) and severely dry (1995) years were investigated on the 
basis of their departure from the mean (Fig. 4a). Results shown in Figures. 4b and c indicate that 
the northern half of the country, which receives more rainfall, suffers more precipitation deficits 
than the southern half. However, Tiffen and Mulele (1994) argue that the impacts of drought are 
felt more in the southern half of the country because of the dominant low precipitation patterns. 
Figure 4b shows that Copperbelt and parts of north-western provinces had the highest 
precipitation deficit during the severely dry year of 1995. During the extremely dry year of 2015, 
the highest deficit was again over Copperbelt and north-western provinces, but it extended further 
to parts of Lusaka and Kariba. In 1992 (extremely dry year), the largest anomalies were observed 
in northern and Luapula provinces over Mpika, Kasama and Mansa. 

It is interesting to note that in both dry and wet years, the rainy season generally started in 
October (Fig. 3) and ended in April or May. These findings indicate that the start and end of the 
rainy season may not necessarily be indicators of whether the year will experience drought; 
forecasting may need to focus on the intra-seasonal variability of rainfall. While the start and 
cessation of the rainy season are fundamental to determining planting dates for the agricultural 
sector, understanding intra-seasonal variability is also fundamental to preparing for probable dry 
spells during the year. The observed consecutive dry days have a huge negative impact on the 
agricultural sector because their occurrence reduces soil moisture (Hottenstein et al., 2015). 
Therefore, studying the variability of consecutive dry days closely will be key for strategic 
planning to ensure food security. 


5 Conclusions 


This study examined the spatio-temporal variability of drought over Zambia. The standardized 
precipitation index (SPI) was used to explore drought occurrence during the period 1960-2016. 
SPI successfully classified extremely dry (SPI value of <—2.0) years as 1992 and 2015, a severely 
dry (-1.9 to —1.5) year as 1995, moderately dry (—1.4 to —1.0) years as 1972, 1980, 1987, 1999 
and 2005, and 26 near normal years (—0.9 to 0.9). The occurrence of consecutive dry days was 
explored and it was found to be strongly negatively correlated with SPI with a correlation 
coefficient of —0.6. 
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Horizontal wind vectors at 850 hPa were investigated during the core of the rainy season 
(December—February, DJF) to ascertain the patterns of flow during years of extreme and severe 
drought. These were contrasted with the patterns of flow in 2007, a generally wet year over 
Zambia. During wet years, Zambia is influenced by large-scale low-pressure zones which boost 
precipitation but tend to be indiscernible during years of dryness. Similar patterns (enhancement 
during wet years and reduction during drought years) were observed for omega (pressure vertical 
velocity) and relative humidity. If followed, these patterns may enable forecasting drought over 
Zambia and aid in strategic planning and management to limit the potential damage caused by 
drought. 
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